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'Hie extent by which the motion of a small particle in a fluid deviates 

from ideal Brownian motion has been the center of considerable interest in the 

past few years. According to the theory of Brownian motion in its simplest 

form*, the velocity autocorrelation function of a small object moving under 

the action of molecular bombardment in a fluid at rest decays exponentially 

'with time. The decay constant is a function of the particle's .mass and radius, 

and of fluid viscosity. This description, which is typical of a Markoff ian 

system with short memory, was believed to be appropriate for Brownian motion 

even for times longer than the velocity relaxation time . Thus it was both 

surprising and exciting when Alder and Wainwright reported a much longer-lived 

velocity correlation function from molecular dynamics computer calculations. 

The long-time persistence of the velocity autocorrelation function, approximately 
- 3/2 

a t decay law, is indicative of previously unsuspected memory effects at 
the molecular level. Subsequent theoretical studies revealed that this long- 
time behavior was consistent with the hydrodynamic description of the particle 
motion. In the hydrodynamic regime, the presistent memory of the velocity 
autocorrelation function has been traced to the existence of local fluid vor- 
tices surrounding the moving object. In fact, even at the molecular level, 
vortex-like patterns of neighboring molecules were observed around test par- 
ticles 5 in computer simulations. This behavior appears to be unrelated to the 

detailed form of the intermolecuiar forces, as more recent computer calculations 

. . 6 

by Levesque and Ashurst demonstrated with a Lennard -Jones interaction potential , 

The first clear experimental verification of the long-time tail of the 
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velocity autocorrelation function was given by Kim and Matta . These 
investigators used a shock wave apparatus to accelerate small latex spheres 


in an atmosphere of argon or air and observed the decaying velocity in a region 
of the shock tube where the superposition of the forward and reflected shock 
waves caused the gas to come to rest locally. The experiment provided conclu- 
sive evidence that the velocity autocorrelation function did not decay expo- 
nentially for all time. The details of the long time decay, however, were 
shown to be rather sensitive to the particle's history prior to the beginning 

t 

of the observation. Thus, the data indicated a long time tail with a decay 

law which was somewhat intermediate between the t" 1 ^ and t ^ power laws, 

- 3/2 

Attempts to exhibit the t ' power law decay by light scattering techniques 
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have yielded less convincing evidence . In these experiments, accuracy is limited 
by the dual requirements that many particles be illuminated and that mutual inter- 
actions between scatterers be minimized. By varying the concentration of scatterers, 
it is possible to satisfy the first requirement. However, it is very difficult to 
obtain a reliable estimate of the effect of particle interactions . 

Over the past year, we have investigated a third approach for measuring 
the non -Markoff ian component in the relaxation mechanism of a Brownian par- 
ticle which combines desirable features of both the shock wave experiment 
and conventional light scattering experimerits . We suggest using the radia- 
tion pressure generated by a G.W. laser to guide an individual spherical par- 
ticle to terminal velocity; at an appropriate time, the beam intensity is 
suddenly lowered to a value at which the radiation pressure is negligible, 
and the ensuing velocity relaxation is measured directly. Recent experiments 
in optical levitation 1 ^ have shown that spherical particles can easily be 
trapped and held stably in a C.W. laser beam operating in the TEM oo mode 
having a Gaussian intensity profile. 

In this report we first review the theoretical background and then 
describe the main -features oi the proposed experiments which arc designed 
to provide information on the velocity relaxation of a small particle. 
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Brownian motion is the apparent If random motion of a small particle 
i: u.terSv’d in a fluid, "Ike oriy/ual invest iga* ions of this pin nemonon were 

made on various objects of colloidal size. The theory, however, has been 
successfully applied to many other situations, from the description of the 
thermal noise of galvanometers, to the motion of ions and polymer molecules 
in solution, even to the influence of rotational diffusion on NMR line widths. 

In its general form, the theory of Brownian motion need not refer to the 
description of a real particle at all, but instead, to the behavior of some 
collective property of a macroscopic system, (e.g., the instantaneous total 
electric dipole moment of a macroscopic sample of fluid. If jt represents a 
typical macroscopic coordinate, its random motion can be described by the 
Langevin equation 

Hfl - -8 i + • <« ' 

where - 3d)f/dt represents the macroscopic damping force and i*(t) is the total 
, external force acting on the appropriate degree of freedom. Usually F(t) is 
composed of various deterministic contributions and of a random part which 
simulates the uncontrollable collision effects. If the random component 
of the force. F(t) lias a very short memory (or if it is s- correlated), the 
velocity correlation function will decay exponentially to zero with a rate 
constant which is given by the damping rate 3: 

^ Y:(t) y ' v(O)>~ e - 0t . (2) 

For a spherical particle in a homogeneous fluid, for example* one finds 

GirnR. 


3 = 


m 


( 3 ) 


where R is the radius of the particle, m is its mass, and n is the fluid 
shear viscosity. 
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A more realistic description of the random motion of a particle in a 
fluid is obtained by the addition of the so-called Basset-Boussinesq non- 
local correction which takes into account drag memory effects. The modi- 

i 

fied Langevin equation takes the form . 1 


-* + t/» 

ni^ = -6irpRv - ^irpR^ - 6R 2 /rnp" j 


dv 

3s fcs 


* r(t) 


( 4 ) 


where p is the fluid density. The first tern on the right hand side is the 
steady state Stokes drag. The second term is the inertial correction asso- 
ciated with the fact that the sphere must displace a volume of fluid equal 
to its own volume before it can proceed forward. The third term is the 
Basset-Boussinesq term which arises from the analysis of the linearized 
Navier-Stokes equation; it introduces memory effects as a result of stream- 
lines induced in the fluid by the particle motion. The last term contains 
possible additional deterministic forces and the usual random force describing 
the collisions of the Brownian particle with the molecules of the fluid. 

If a particle is injected in the fluid with a non- zero initial velocity, 
the velocity relaxation to equilibrium is described on the average by an 
equation with the same formal structure as the modified Langevin equation (4) 
without the random force (since- thdiaverage random force is zero). The linearity 
of this equation also insures that the velocity autocorrelation function will . 
relax to its equilibrium value according to the same equation as that satis- 
fied by the average velocity. The key point in this discussion is that the 
presence of the non- local term-causes a significant departure for long times 
from the exponential decay which is expected for both the average velocity 
autocorrelation function according to the simple Langevin description. Thus , 
we can obtain the long time behavior of the autocorrelation function by measuring 
the long time behavior of the average velocity, from Ik}. (4) we see 
that if < 0 for all s < t, the Basset term causes a delay in the damping 
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process so that the particle drift velocity approaches aero over a Longer 
time scale than the usual dtaractensuc damping time w U^nR) of ordinary | 

Brownian motion. 

This is not very surprising from a hydrodynamic point of view (the 
streamlines of the perturbed fluid curl around the moving sphere and "push" : : 
it from behind in the direction of the motion) . What is surprising is that 
Eq. (4), in spite of its complicated structure, can be solved exactly if the 
sphere’s initial velocity is constant, v(t < 0) = v Q , and if the mean value 
of the external force is also constant, <F(t))> = F. First we consider the 
case F = 0. Because Eq. (4) is linear in the velocity v, the average drift j 
velocity of the Brownian particle satisfies the same equation as (4) without 
the random driving force. In the following development we will have no further 
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need to consider the random velocity v so that without danger of confusion, the 
symbol v(t) from now on will denote the average drift velocity. Its equation 
of motion is obtained by taking the average value of each term in Eq. (4): 


where 


and 


dv 

at 


A i 

♦ Pv +YJ ds 


. _ +6ffnR 
p ” m+m^ v ; ■ 

v = bR 2 (Trnp) 1/2 
' m+m£ 


m _ 2 tt n 3 
Ifif — j R : p : 


X L = 0 
ds 


(5) 

( 6 ) 

(7) 

( 8 ) 


On the surface ft might appear that the drift motion of the Brownian particle 
is controlled independently by its Stokes damping rate 3 and by the strength 
y of the memory term. In reality, the memory effect in the context of this 
model imposes a kind of "universal" deviation from the usual Brownian motion 




if the time scale of the evolution is measured in units of fi' 1 . We can see, 

. ' ’ • ' ■ ' . . . . . . . : 

this clearly if we introduce the dimensionless time S <*\et in Eq. (5). The 
evolution equation becomes 


dv 


a? 


L- + v + 



€ 

dV 


dv i 

dT' /c 2 





Hence the extent' to which the Basset .term affects the particle Velocity 
controlled uniquely by the parameter f 


(9) 

is 


JL- Jl 

/p nr 

where p is the mass density of the sphere. Surprisingly enough, the scaled 

if -' - ■ • - 

equation is no longer governed directly by the fluid shear viscosity n so that 
different fluid media with the same mass density p f will have the same effect 

i r 

on the particle velocity. This is, of course, true only if the observed time 
dependence of the drift velocity is analyzed on a time scale such that the 
unit of time corresponds to one Stokes decay time. 

Equations (9) and (10) also indicate that in order to enhance the devia- 
tion from ordinary Brownian motion, one should use a fluid medium with a 
relatively high mass density. This, in fact, is not entirely feasible because 
high fluid density also implies large shear viscosity and short Stokes relaxa- 
tion times . .Thus , time resolution problems may become serious if one should 
consider a fluid-like water where the Stokes relaxation time for 10 y m radius 
particle is only of the order of 20 ysec. 

By contrast, a 10 pm radius particle in air at a pressure of 1 atmosphere 
and a temperature of 20°C relaxes to equilibrium in a time of the order of 1 
msec*, which is easily measured, e.g,, by laser Doppler spectroscopy techni- 
ques. However, as noted above, tii^ ; -fi^tivo''-importande of memory effects 
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decreases when the fluid density is small. Thus , in the case of 10 yn glass 
heads in air at standard conditions, one finds » 4 x 10"“, 

We now consider the solution of Eq. (9) in order to put our comments on 
a more quantitative basis. The exact solution of Eq. (9) can be derived by 
a rather straightforward application of Laplace transform techniques, The 
velocity decay is rather complicated: 


v(g) _ y+/ii Z “l 


v_ ZV\l 

\ y /y^-T 


exp [(p-.'5^T)2 {] 
exp [ Cy+4i T -i) 2 5] 


2^-1 


where I 

i / 




erfc [(y-/y z -l)/£ ] 
erfc [(y+/y 2 -l)/f ] 


( 11 ) 


( 12 ) 


We can obtain the long-time behavior of the Velocity by performing an asymp 


totic analysis of Eq. (11) in the limit § ■+ ® (since 5 = 3t) with the result 

vCg)^ r (15) 


We are thus able to demonstrate mathematically how collision memory effects 
lead to an inverse power law decay in the velocity relaxation over long times, 
which should be compared to the exponential decay of simple Brownian motion. 

Up to this point, we have considered decay in the absence of external 
forces. However, if the external forces F in Eq.(4) are constant, we can 
easily obtain the velocity decay from Eq. (11) by rescaling the velocity. 

In fact, we simply subtract the quantity T = F/[(m + ni^) 0 2 ] from v(g) and v Q . 

In view of the size of the memory term parameter we have also 

developed an accurate perturbative solution to first order in y//$ , which 
although not limited to constant F, can be applied to the case F = 0. After 
taking a Laplace transform of Eq. (5) we find 
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v(s) 


1 +’ y V v^T 
s+X+y'/^ 


(14) 


where 


Hr * 


7(5) » L(v(0) 

To first order in y' , Eq, (14) takes the fonn 


7(s) 


+ °Cy’ 2 ) 

i, Eq. (17) can be easily inverted with the result 

vT 

>-= e / 

/A I 
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•o 


SM-e- 5 e-« 

V - /P 


( 16 ) 


(17) 


J % A 2 dX - &5 1/2 ^(2,3/25-0 


(18) 


The second term on the right hand side of Eq. (18) is the so-called Dawson 
11 

function ; the third is the better known confluent hypergeometric function. 

It is clear that for short times (£ <_ 1) the dominant contribution is the 
exponential term that describes the usual relaxation of a Brownian particle 
(this is especially true because of the smallness of the term y//$ )V For 
long times instead, a careful asymptotic expansion of the Dawson function 
and of the confluent hypergeometric function leads to the following limiting 
behavior of v(?) • 

1/2 d + ?" l + ••*) f • (19) 

Clearly , the asymptotic behavior of the drift "velocity for a particle subject 
to a zero external force (for example, in a gravity -free environment) is dras- 
tically different from an exponential. There remains to see at what stage of 
the evolution a clear depature from the memory-free behavior becomes observable. 
This is shown in Fig. 1 for several values of the ratio y//3 . 

The experimental tests which we have proposed to study memory effects in 
an earthbounu jajuuatory, envision, the suspension of a small glass bead in a 
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powerful (a few watts) focused C.W. laser beam. 'Hie particle is initially at 
rest in air at various pressures and can be trapped stably in the restoring' 
radiation pressure field provided by the Gaussian laser beam. If the beam 
intensity is suddenly changed by a small amount, the particle will seek a 
new equilibrium position undergoing a modification of the usual damped har- 
monic motion. Hie modification will come about as a result of the memory 
term in the equation of motion. 

Hie Langevin equation for the drift velocity of a particle subject to 
viscous damping and a harmonic restoring force lias the form 



where K is the usual dimensionless time, k is the scaled restoring force 
constant and z is the instantaneous distance of the particle from equili- 
brium. Unlike the case of a free particle where an exact solution of the 
Langevin equation can be obtained without much effort in this case, the 
Laplace inversion requires the solution of a cubic equation. It is clearly 
preferable to forgo the search for an analytic solution, and to take advan- 
tage instead of the smallness of y//$ to obtain a perturbative solution 
which is correct up to terms of order (y ) . We consider again the 
Laplace transform of Eq, (20) and eliminate the variable z (?) with the help 
of -the relations 



s z(s) - z D «.!v(s) 

where z(s) 5 L (z(c)) and z = z(§ - 0). After some minor algebraic steps 
we find 


io 


vis) 


k zj s 



5 ’'o' 11 Z 0 
s 2 + s +k 


+ Y -S 


[k + v ) s + k y, 
t/T 1 0 0 0 


( 22 ) 


(s 2 + s + k) 2 


where y ! 


v/iT 

?f * 


'Hie inversion of Eq. (22) is rather laborious, but 


straightforward. The result is 
v(£) - A e"^ + Be 


2y C 

/r 


r r 

.P e " P5 J e X dX - ^ o' 5 J 

0 0 


/QC 


e x * d\ 


+ ^ ? 1/2 (D ^ (2,2/2* -PC) + E ^ (2,3/2; - QS)] 
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Before we proceed to identify die numerous parameters of this solution, We 
must point out at once that the last two terms on the right liand side of 
Eq. (23) represent the memory correction. Furthermore, the evolution of 
v(C) proceeds according to two distinct time scales set by the rates P and 
Q. We now discuss the various symbols appearing in Eq. (23) . The parameters 
P and Q are solutions of the quadratic equation s 2 + s + k = 0 (see Eq. (22)), 


x.e. , 


P = J (1 + /l - "4 k ) 
1' 


(24) 

(25) 


Q = • J (1 - /l^4k ) 

They can be real or complex as with the ordinary problem of a damped harmonic 
oscillator. The constants A and B which govern the ordinary (non-memory 
related) part of the solution are given by 


A = 


B = 


kz + P v 
o o 


TV 

1PF 


k z o + Q v o 


(26) 

(27) 
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!! 


Note that In the typical experimental situation which is being considered we 


iwn’t to have v ” 0 
* e 


h z + v - 2k v 
c n - 1° 0 0 


CP - Q ) 3 


k v, 
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E s 
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z rt + vj 

0 o' 


CP - Q ) 3 

k v 0 - Q(k z Q * v 0 ) 


CP - Q ) 2 


C28) 

(29) 

(30) 


We can easily verify that if k 0, Eq. (23) properly reduces to Eq. (18) 

t 

(force free case), Hie long time limit of Eq, (23) however provides us with 

a surprising result which helps in dispelling some of the early attempts in 

* - * / 2 

the literature to associate excessive fundamental significance to the £ 

behavior of the long time force free ease. In fact, the long time limit of 

• * 

Eq. (1.6) is seen to be^of the form 


v(0 ~ 


1 1 
Q " p 





(31) 


The power law is now of the type s'* 5 ' 2 indicating that the asymptotic behavior 
is apparently a sensitive function of the external field of force. 

We should note here that while the solution of Eq, (23) converges to the 
one obtained in the force free case for k ^ 0, it is no longer possible to * 
try to compare Eq. (31) with its long time counterpart in the k + 0 limit. 

This is because as k 0, the parameter Q approaches zero thus violating the 
asymptotic conditions for Pf; and which are necessary in. order to derive 
Eq. (31). 

It is interesting to observe that in the overdamped case (k < 1/4) , three 
phases of the relaxation should be investigated carefully: 
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a) r$*miQ5«l 

In tlus case, the motion is essentially indistinguishable 
from that of an overdamped oscillator 

» 

b) P$ > 1 but Qs < 1 

Here, one should be able to detect evidence of memory 
effects 

c) P5 » 1 # Q5 » 1 

Ihis corresponds to the as>Tnptotic decay region which is 
governed by the power law (Eq. (31)). This situation appears 
difficult to study experimentally because of the long time 
scale involved. 

In addition to the exact and perturbative solutions to special cases of 
the modified Langevin equation, we also constructed a numerical solution which, 
is applicable to arbitrary forces F (t) . This solution was based on the Adams - 
Bashforth-Moulton predictor/correetor technique for solving differential equa- 
tions. In extensive tests between the numerical and analytical solutions, 
agreement to four significant digits was demonstrated throughout thq region 
of interest, The numerical solution can be used to predict the motion of a 
levitated glass bead when the laser undergoes small sinusoidal fluctuations 
in power. By varying the frequency of these fluctuations , the resonance of 
the system can be probed to yield critical information on the size of the 
memory term. This point, however, needs additional study, and will be 
analyzed in subsequent progress reports . 

We discuss now a few possible ways to verify the predictions detailed 
above. We are especially interested in measuring the velocity relaxation 
of a small glass sphere (diameter of the order of 10-40 ym) in a gas. In 
an earthbound laboratory, a possible way to accomplish this goal is to 
levitate the particle using the radiation pressure produced by a focused 
Argon-Ion laser having a Gaussian intensity profile. As shown schematically 




in Fig. 2, the laser ‘"am is deflected upward and focused on a collection 
of particles resting on a microscope slide, A cylindrical piezoelectric 

■ ■ . ' ' A\> 

crystal is used to break the Van der Waals attraction between the glass 
beads and the slide. With some careful handling, a sphere can be trapped 
in the laser beam and propelled upward above the focal plane of the lens 
where the effects of gravity, radiation pressure and radiometric forces (if 
present) balance each other out. The beam profile is 'such that the suspended 
spheres are trapped stably as a result of radial restoring forces if sufficient 
precaution is exercised in shielding the system from external disturbances 
(air drafts, mainly). Once the particle is levitated, the laser power is to 
be suddenly increased or decreased to a new level by means of an electro -optic 
device and the ensuing motion will be measured, 

A direct recording of the motion can be done with a high speed streak 
camera of the rotating drum type. The required resolution of 0.1 ysec appears 
to be available commercially. At a subsequent time, we will develop an optical 
heterodyne detection scheme for measuring the particle velocity directly 
through the Doppler shift of the backscattered light. Real-time frequency 
measurements of the beat signal between a reference beam and the Doppler 
shifted one will give direct information on the relaxation process of 
interest. These measurements appear quite feasible because of the high 
intensity of the backscattered radiation, but it is difficult to predict 
how accurate the readings can be for the required, long times . Of special 
interest is the analysis of the scaling relation developed in the theoretical 
section* For this purpose, we will modify the cell design (item 6, Fig. 2) 
so that the atmospheric pressure can be varied in a controlled way. This is, 
of course, quite easy to do because the required operating pressures will be 


of torr). Under low , »’essure conditions, the particle motion towarj the 
perturbed equilibrium coni ^.urat ion should be underdamped and easy to follow 
with moderate time resolution. 

It appears interesting to force a sinusoidal modulation of the laser y 
intensity at various frequencies and to observe the particle response. As 
the modulation frequency is varied across the particle resonance, additional 
information on the non-Markoff ian components of the velocity will be sought. 

Certain aspects of the velocity relaxation cannot be measured in an 

r 

earthbound laboratory owing to the presence of gravity. In particular, the 
t power law derived previously in this paper applies only to velocity 
relaxation in the absence of forces. It is this relaxation which we will 
measure in a gravity-free enivronment . Now, the laser can be used to propel 
a sphere to terminal velocity in the scattering cell. After the laser is 
reduced to the level of a probe signal, the velocity relaxation will be 
measured using the same heterodyne detection system. Other systems for 
zero-gravity experiments may also be developed. Since there is no gravita- 
tional force to overcome, low power lasers may be used to propel the spheres. 
In addition, the spheres may be trapped by configuring two laser beams of 
equal intensity antiparallel to each other. 


log ( t(£)/v(0) ) 
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